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Abstract 

We consider Yang-Baxter equations arising from its associative analog and study cor¬ 
responding exchange relations. They generate finite-dimensional quantum algebras which 
have form of coupled GL(A) Sklyanin elliptic algebras. Then we proceed to a natural gener¬ 
alization of the Baxter-Belavin quantum 77-matrix to the case Mat (A, C)®^(8)Mat(M, C)®^. 
It can be viewed as symmetric form of GL(AM) 77-matrix in the sense that the Planck 
constant and the spectral parameter enter (almost) symmetrically. Such type (symmetric) 
77-matrices are also shown to satisfy the Yang-Baxter like quadratic and cubic equations. 
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1 Introduction 

The associative Yang-Baxter equation 


a 

a 

a 

a 


0 

[3 


R12R23 — R13R12 + R23R13 


( 1 . 1 ) 


appeared in P as coassociativity condition (A (g) id^i) o A = (id^i (g) A) o A for the (principle) 
derivation A ■. A ^ A® A defined in an associative algebra A by A (a) = aR — Ra, Ya & A and 
some R & A® A. In this paper A = Mat (A, C), i.e. N x N matrices over C. 

The equation fll.ll) was extended for R depending on additional parameters [H] and used for 
description of structures behind the classical Yang-Baxter equation on elliptic curves and their 
degenerations [6]. Let us write it as 

RtRl = RI,RY + R®Rl, R%,^RL{^a-‘t), (i-2) 


(see (1A.23P ) where Za, Zb, Zc, h,ri denote free generic complex parameters. Indices a,b,c are 
distinct numbers of tensor components in Mat(A, C)®”-^. is dehned in a-th and 6-th com¬ 
ponents. If c 7^ a, 6 then acts on the c-th component of Mat(A, C)®*^ by A x A identity 
matrix. 


The elliptic solution of fll.2p was found in 
quantum i?-matrix mu (in vector representation), 
quantum Yang-Baxter equatioiii]: 


It is the Baxter-Belavin’s elliptic GLat 
T he latter means that satishes the 


JDfi JdR _ JdR JDfi pft 

^ab^ac^bc — ^bc^ac^ab 


(1.3) 


The property (11.211 of the Baxter-Belavin i?-matrix was rediscovered in our papers [IHl [14] in the 
framework of integrable systems and related topics - Painleve equations, Schlesinger systems 
and their quantization via the Knizhnik-Zamolodchikov-Bernard (KZB) equations. 

^See [To] for reviews on the Yang-Baxter equation and related structures. 
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Elliptic function identities. The associative Yang-Baxter equation fll.2p in the scalar case 
(when N = 1 and a,b,c = 1, 2, 3) is the Fay trisecant identity on elliptic curve 

4 >{zi-Z2, h) 4 ){z 2 -Z3,r]) = (j){zi-Z3,r])(j){zi-Z2, h-T]) + (j){z2-Z3,r]-h)(j){zi-Z3,h) (1.4) 


The function (j){z, u) satisfying this equation is the Kronecker function. It is dehne as 




'&'{Qi)'d{z + u) 
'd{z)d{u) 


(1.5) 


where d{z) is the odd Riemann theta-functiorj^ flA.3p . In this respect the Baxter-Belavin’s 
R-matrix can be viewed as non-abelian generalization of the Kronecker function. In [IT] we de¬ 
scribed a list of elliptic function identities and properties together with their R-matrix analogues. 
In particular, it is easy to verify that the Baxter-Belavin i?-matrix dehned a^ 


+ h)Ta® T_„ , Res ^ 12 ( 2 :) = NP 12 , 
' 2=0 




see 


dehnition of (^^(z, (Uq, -t- h) in (IA.22P ) satishes the skew-symmetry property 


TDfl _ C> — 

^ab — ~^ba 


and the unitarity condition (see Appendix) 

RlRl = 1 ® 1 N^MNK) - p(z.-2,)), 


( 1 . 6 ) 


(1.7) 


( 1 . 8 ) 


where p{z) is the Weierstrass p-function. Condition fll.7p is analogue of (j){z,u) = —(j){—z, —u) 
while fll.Sp is similar to flA.7D (j){z^u)(f>{z, —u) = p(z) — p(u). 

In trigonometric (hyperbolic) or rational cases the Weierstrass p-function entering fll.81) 
equals p(z) = 1/ sinh^( 2 ;) or p(z) = 1 /z^ respectively. The simplest example of the R-matrix 
satisfying fll.2p . fll.3p and fll.7p . fll.8p is the quantum Yang’s R-matrix [23]: 



1 (g) 1 YR 12 

h z 


(1.9) 


It is the R-matrix analogue of function l/h-|- 1/;^. The class of R-matrices under consideration 
includes nontrivial trigonometric and rational degenerations (of fll.bp h which satisfy fll.2l) . fll.7p 
and fll.8p . These type R-matrices can be found in [2] and [22l[T5] for trigonometric and rational 
cases respectively. In GL 2 case the corresponding 7-vertex and 11-vertex R-matrices were found 
in [7j. See also [B] which results are (presumably) gauge equivalent to those of [21 1221IT3] . 


Purpose of the paper. 


1. Yang-Baxter equation with two Planck constants. As we will see it follows from fll.2p . 
fll.7p and fll.8p that the quantum R-matrix satishes also cubic equations of Yang-Baxter type. 
The one closest to the original one fll.3l) have the following form: 

+ RizRuRu . ( 1 - 10 ) 

^In trigonometric (hyperbolic) or rational cases (j){z,u) = coth( 2 ;)-|-coth(u) or (j){z, u) = I/z-I-I/m respectively. 
•^In [T3] this i?-matrix was found in different - Richey-Tracy form [12] . 
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where again = R^bi^a — ^fe)- When h = rj fll.lOp coincides with the Yang-Baxter eqna- 
tion fll.dp . Another (more general) type of eqnation contains p-fnnction entering the nnitarity 
condition fll.Sp : 


“ R’kRaRvz = Ra"MNn) - p(Nh)) . (1.11) 

Let ns mention that eqnations of (II.lip type were considered in [20] in the context of donble 
Lie (and the qnadratic Poisson) strnctnres. 

2. Quantum algebras. Given a qnantnm i?-matrix one can dehne the following (Sklyanin 
type) algebra nsing exchange relations: 

RUz - w)L^,iz, S)L^,{w, S) = L^(n;, 5)L^(^, S)R’l,iz - w) (1.12) 

with L^{z,S) = tr 2 j • For N = 2, i.e. in the Baxter’s case eqnation fll.l2p provides 

the Sklyanin algebra IZH for the components Sa, S = ^^Ta. In particnlar case (see (13.Op ) 

a 

snch L-operator is jnst the i?-matrix (II .Op itself, and the exchange relations (I1.12p are fnlhlled 
identically dne to Yang-Baxter eqnation fll.3p . It means that fll.l2p follows from fll.3p by treating 
one of tensor components as some Sklyanin algebra modnle. 

In the same way we can dehne exchange relations corresponding to fll.lOp as follows: 
Rl,{z-w)^{z, S’^) + RJl,{z-w)L’l{z, S’^)Ll{w, S^) = 

(1.13) 

= lI{w,S'')L’1{z,S’<)rUz-w) + LV,(w,S’>)L'S{z,S’‘)IiUz-w) 

Again, when h = rj it coincides with the ordinary exchange relations fll.121) . As we will see this 
eqnation dehnes conpled Sklyanin algebras, i.e. given a pair of Sklyanin algebras for and S'^ 
the eqnation (11.131) provides commntation relations between and S^. 

3. Symmetric i?-matrix in Mat(Y, C)®^ ® Mat(M, C)®^. Here we consider the following 
extension of the Baxter-Belavin i?-matrix (11.61) : 

5 ^ 12 , 12 (^ 1 ^)= ^ exp{ 2 mN^h)R^ 2 iz + (I.I 4 ) 


where Tq is Tq of size M x M. We call this i?-matrix symmetric in the sense that both argnments 
(the Planck constant h and the spectral parameter z) are averaged over lattices Zat x Ztv and 
Zm X Zm respectively. In the simplest rational case corresponding to Yang’s i?-matrix fll.Qp the 
expression fll.ldp is rednced to 


5^12,12(^1 h) — iff 


In G Iv <8) P12 , P12 G 1m G) 1m 


h 


+ N- 


(1.15) 


The limiting cases Y = 1 or M = 1 are given as follows: 

5^12,12(^1 

M=l^ 

R&(z) 
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(1.16) 






























We will show that such an i?-matrix satishes a set of relations similar to those for the ordinary 
/^-matrix. In particular, we have the following generalization of the associative Yang-Baxter 
equation fll.4p : 


^ 12,12 ^23,32 


~ ^13,32 ^12,13 


^23,31 ^13,12 ) 


(1.17) 


where ~ ^abdii^a, — Zb, ha — h^). The cubic equations arise similarly. For example, the 

generalizations of fll.llj) takes the form: 


(1.18) 


^ 12,32 ^ 13,13 ^ 23,32 ~ ^ 23,13 ^ 13,32 ^ 12,13 + ^ 13,12 yPi^h^2) Pi^hl^)j , 

^.32,12^13,13^32,23 ~ ^13,23^32,13^13,12 + N M 91]^2,i3 (^^("^^ 23 ) ~ ^(^ 2 : 13 )^ . 
where Zab = Za — Zb and = ha — hi- The Yang-Baxter like equation fll.lUp is generalized as 


follows: 

^ 12 , 32 (^ 12 ) ^ 2 ) ^ 13 , 13 (^ 13 ) ^ 13 ) ^ 23 , 32 (^ 23 ; ^ 32 ) + ^ 12 , 32 (^ 12 ) ^ 13 ) ^ 13 , 13 (^ 13 ) ^ 32 ) ^ 23 , 32 (^ 23 , ^ 13 ) 
~ ^23,13(^23, ^13)^13,32(^13) ^32)^12,13(^12) ^13) + ^23,13(^23) ^32)^13,32(^13) ^13)^12,13(^12) ^32) 


and 

^32,12(^32) ^12) ^13,13(^13) ^13) ^32,23(^32) ^23) + ^32,12(^13) ^12) ^32,13(^13) ^13) ^32,23(^13) ^23) 
~ ^13,23(^13) ^23)^32,13(^32) ^13)^13,12(^13) ^12) + ^13,23(^32) ^3)5^32,13(^13) ^13)^13,12(^32) ^12) 
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15-51-52031 HHCq. The work of A.L. was partially supported by Department of Mathematics NRU 
HSE, the subsidy granted to the HSE by the Government of the Russian Federation for the implemen¬ 
tation of the Global Gompetitiveness Program and by the Simons Foundation. The work of A.Z. was 
also partially supported by the D. Zimin’s fund ’’Dynasty”. 


2 Yang-Baxter equations with two Planck constants 


Consider a quantum i7-matrix which satisfies the associative Yang-Baxter equation (11.21) and 
the properties fll.7p . fll.Sp . Then it is easy to get a set of cubic (in R) relations. 


Proposition 2.1 

(El; 


The following cubic relations are valid for a common solution of lil.S\} . |fi. ?[), 


pb pb pb _ pb pb pb 11 

-'l'12-ll'13-'^23 ~ -ll'23-^13-^12 ) 

RlRlRl + RlR'‘,Rl = -N^f.'(Nh) 1 . ® U 0 1 ,, ( 2 . 2 ) 

RlRlK - RLRl,Rl = Rii" MNri) - p(nk)) , ( 2 . 3 ) 

KiRlK + RlRLRl = RlRlK, + RlRLRl . ( 2 . 4 ) 


where = Rab{za — Zb) and a, b, c are distinct numbers from the set {1, 2, 3}. 
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Proof: First, notice that the primary one is the third identity fl2.3p since all others follow from 
it. Indeed, the Yang-Baxter equation fl2.ip follows from fl2.3p in the case h = rj. 


The second equation (I2.2p appears from (I2.3p in the limit rj —)■ —h: one should use the 
classical limit flA.25p . and then the skew-symmetry property (11.7p to arrange the indices in a 
cyclic order. This gives fl2.2p up to permutation of indices. 


The fourth relation (12.4p can be called Yang-Baxter equation with two Planck constants 
because its structure is similar to fl2.ip and it is coincide with fl2.ip for h = rj. It is easy to see 
that fl2.4p follows from fl2.3p as skew symmetry of fl2.3p l.h.s. with respect to interchanging h 
and rj. 


Thus we need to prove 
side: 


Consider fll.2p and multiply both parts by ^ from the left 

(2.5) 


]D^~V jdV _ _L 

■^6c ^ah^hr ^hr ^ar.^nh ' ■^h'^ -H., 


^ab-^ ^bc 


^bc 


^ab 


^bc 


^bc 


Now interchange the indices 2,3 in fll.2p and multiply both parts by R^^ 


r>h jdV jdV _ jdV D^—>7 JdV _ 

^ac^cb^bc ^ab^ac ^bc ^bc ^ab^bc ’ 


( 2 . 6 ) 


where in the last term we have already used R\ ^ = —R^^ From fll.7p . fll.Sp it follows 
that and Rf^^Rfl^ are scalar operators. Subtracting (12.6p from (12. 5 p we get (12. 3 p with 

h ■.= h — rj. ■ 

Let us also remark that the identity fl2.2p appeared to be related to i?-matrix valued Lax 
operator for the classical Calogero-Moser model [13]. Consider the following block-matrix: 


n 

Sab)Rabi^a “ Zb) , (2.7) 

a,6=1 


where Eab is the standard basis in Mat(n, C), i.e. {Eab)cd = ^ac^bd- Then the diagonal blocks of 
trC^lh) are scalar operators, and the scalar functions are obtained as if C{h) were element of 
Mat(n, C) with matrix elements Lb{h) = (1 - 6ab)N(f){Nh, Za - Zb): 


Y(fi)).„ = i®...®iO‘(fi))„„. 


( 2 . 8 ) 


When n = 2 this equation is equivalent to the unitarity condition fll.Sp . while for n = 3 fl2.8p 
reproduces (12.2p . 


One more application of fl2.2p comes from the classical limit flA.25p . The identities which 
appear from (12.2p together with (IA.29P provide sufficient conditions for compatibility of the 
KZB connections (see details in [T5]L 


Finally we conclude that the identity (12.3p is of great importance. On one hand it reproduces 
the Yang-Baxter equations, and in this sense it ’’knows” about quantum integrability and related 
algebraic structures including quantum groups, Sklyanin algebras, e.t.c. On the other hand 
(12. 3 p also ’’knows” about classical integrable system of Calogero type. At last, the same identity 
provides the quantization of the Schlesinger systems via compatibility of the KZB connections. 
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3 Quantum algebras 


In this section we discuss differen types of quadratic finite-dimensional quantum algebras arising 
from exchange like relations. The quantum (Lax) L-operator [211 IE] is defined as 


L\z,S) = L\z) = ti 2 {RU^)S 2 ), 5 = 

a 


(3.1) 


where {Sa, cn G is the set of generators of the quantum algebra A. For the elliptic i?-matrix 
(II-bh we have 

= '^T^Sa(Pa{z,UJa + • 


For the vector representation of A given by 

p^{S^) = T_^ e Mat(A^,C) 


(3.2) 


(3.3) 


the Lax operator fl3.2p coincides with the i?-matrix fll.61) . 

Finite Heisenberg group (see flA.lSD in Appendix) is the simplest example of quantum algebra 
which comes from i?-matrix relations for the L-operator fl3.ip . 


Proposition 3.1 Relations 

L’l{z)^{w) = Ll-^^{w)R^,,{z -w)- R-,\z - w)lI^^{z) (3.4) 


or 

^{w)L’l{z) = R^ 2 iz - - l1'^'^{z)R'^^{z - w) 

with L-operator defined in iSA) and R-matrix U.6\) are equivalent to 


S, S 2 = NP,2 S, , 


(3.5) 

(3.6) 


i.e. 


S^Ss = 


Kn 


a,p 


(3.7) 


in components. The associativity condition for the triple product L\{z)VI{w)L\{x) follows from 
the associative Yang-Baxter equation di.il)) . 


The proof of fl3.4l) . fl3.5p follows directly from the Fay identity fll.4|) . The statement about 
associativity is analogues to the derivation of fll.ip in [T]. 

Notice that the commutation relations 

|if (j), LUw)] = lLl*”(z), r:;(z - «,)! + |i5+’(u>), rU^ - »)l (s.s) 

are equivalent to those for the Lie algebra g\^: 

[S,,S2]=N[P,2,S,]. (3.9) 

Obviously the same expression (r.h.s. of fl3.8p l can be used for definition of the Poisson-Lie 
brackets {S'!, S 2 } = A^[Pi 2 , Si] on gl^: 

{l;(2), lkw )} = il;+''( 2), ii3'(z -!»)] + iiy’'(m), ru ^- - w)]. ( 3 . 10 ) 
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with the classical matrix-valued function L{z, S) on the phase space gl^ parameterized by Sa 
{S = ^aTa). The Jacobi identity is due to (II.2^ . The relation (13.1 Oh differs from the custom 
classical exchange relations 


{li{z),l 2 {w)} = [h{z),ri 2 {z - w)] [l 2 {u!),ri 2 {z - w )], 


(3.11) 


where the classical r-matrix ri 2 {z) (1A.25P is used, and the Jacobi identity is fulhlled due the 
classical Yang-Baxter equation flA.26p for ri 2 {z). Alternatively, we can say that in the case 
under consideration the constant h (and rj) is not the Planck constant entering the quantum 
i?-matrix but rather additional spectral parameter entering the classical r-matrix. Indeed, the 
hnal commutation relations (13.9h are independent of h, rj (as well as they are independent of 2 ; 
and w). Such interpretation is close to the consideration suggested in [T8] . 

Exchange relations with two Planck constants and coupled Sklyanin algebras. Con¬ 
sider a pair of the quantum Lax operators 


and 


i*(;) = s«) = + h) 


L^z) = L^iz, S'')=J2 tJ'Mz, + r,) 


(3.12) 


(3.13) 


Each of them dehnes its own Sklyanin algebra via exchange relations (11.121) . See details in 
the Appendix. Here we suggest another exchange type relation which provides commutation 
relations between and S^. It is of the form: 


R’l2{z-w)L'l{z)L2{w) + BJl2{z-w)L^{z)Ll{w) = 

= L^2MLUz)RUz-w) + L:>{w)L^,{z)RUz-w) 


(3.14) 


Proposition 3.2 Equation Ii3.14\ ) for the Lax operators L^{S^) and L^{S^) is equivalent to the 
following commutation relations for the components of and S^: 


ii-,+ h'.7 =0 • 


where the structure constants f^’^ ^ are given by 

for /3 7 ^ 0 : fj^ .^ = Ei{u^ + h) - Ei{Ua-p-^ + V) + + V) - ^l(^/3+7 + ^) , 

for /3 = 0 : f^’J] = p{uj^ + h) - piuja-j + p) ■ 

Proof: Here we use short notation p>^{z) = ip^^z^oja + h). Compute 

RUz-w)LI{z)L^2{w) - LI{w)L\{z)R:12{z-w) =Y,Ta®TpX 

a,0,'y 


(3.15) 


(3.16) 


















By redefining the snmmation index for the second term asy—)-a —/9 —ywe get for the following 
coefficient behind the tensor component 

«>/3,7 


= T„ ® 4?,-, . 

«,/3,7 

where we assnme (po{w) = 1. The last eqnality follows from flA.Sp . flA.Qp . In this way we 
reprodnce the hrst term in fl3.15|) . The second term in fl3.15|) is obtained in the same way by 
interchanging h and r]. ■ 

Notice that when h = rj the structure constants f^’^^ = coincide with those for the 
Sklyanin algebra flA.3ip . 


Recently the coupled Sklyanin algebras appeared also in different way - via the modular 
double [^. In that case two algebras depend on different modular parameters r. In our case 
this parameter is the same but the Planck constants are different. 


Remark. One can also consider (quasi)classical limits of the Yang-Baxter 02.41) and exchange 
relations 03.141) . There are different possibilities for the limits since we deal with two Planck 
constants. In particular, when h —)■ 0 (while r] is hnite) 02.41) is reduced to 


[Ri2, R 23 ] + [-^ 13 ^ 23 ] + [^"12, -R 13 ] — 0 , 


(3.17) 


where Vab is the classical r-matrix 0A.25D . At this stage the generators become classical 
variables Sa while are still quantum. Then we get a ’’half-quantum” Poisson structure 
between the commutative classical variables Sa and noncommutative variables S^. Taking then 
limit r] ^ 0 the equation 03.171) is reduced to the ordinary classical Yang-Baxter equation 
0A.26I) . and the Poisson structure is the standard classical Sklyanin algebra. We will describe 
the quasiclassical limits in our next paper. 


4 Symmetric i?-matrix in Mat(A7, (8) Mat(M, 


In this section we consider an extension of the i?-matrix 01.61) in Mat(Y, C)®^ to the matrix 
function in Mat(Y, C)®^ 0 Mat(M, C)®^. This type of R-matrices (see e.g. [TTl [161 E]) appears 
naturally when N := NxM from the original Mat(Y, C)®^-valued one. Our definition is slightly 
modified in a way to get expression which is symmetric in spectral parameter z and the Planck 
constant h. The simplest example of the i?-matrix which we are going to discuss here can be 
given in the rational case. Consider the following expression: 


^12,12(^1 Z 2 ,hi h ^) 


^ l^r 0 Iat 0 P12 

h - h 


_l_ P 12 0 1 m 0 1 m 

2^1 - ^2 


5^12.12(0 - ^2, - h) e Mat(Y, C )®2 0 Mat(M, C)®^ , 


(4.1) 


where P 12 is the permutation operator in Mat(M, C)®^, and l^v (or 1m) is the identity matrix in 
Mat(Y, C) (or in Mat(M, C)). 93^2 12 tensor indices. The hrst pair is for the components 
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of Mat(A^, C)®^, and the second pair is for the components of Mat(M, C)®^. When M = 1 fl4.1|) 
coincides with the Yang’s i?-matrix (II .hh ~ ^ 2 ) £ Mat(Y, C)®^. In the same way when 

N = 1 gH) gives - h) e Mat(M, C)®^. 

The expression (14. Ih is not a new i?-matrix of course. Multiplying it by Iat 0 Iat <8 ) P 12 we 
get a special form of the Yang’s i?-matrix in Mat(YM, C)®^. 

Below we define elliptic analogue of fl4.ip . But hrst we recall some important properties of 
the Baxter-Belavin’s i?-matrix which will be to extended to the case Mat(Y, C)®^0Mat(M, C)®^ 
for the proper definition of underlying set of functions. To be exact, our goal is to define a set 
of functions ^a,a{z, h), a ^ Zn X Zn, a ^ 'Em x Em which are invariant with respect to shift 
by the full lattice periods a^a + N,a^a + Nr and a ^ a + M, d —)■ d + Mr. 

For the reasons given above (and for the brevity sake) we will refer to 94^2 12 tyP® expression 
as symmetric i?-matrix though it is not such symmetric in the elliptic case as in fl4.ip . 

Symmetries of Baxter-Belavin i?-matrix. The quantum elliptic i?-matrix (II.6 p is Zat x Zat 
symmetric, i.e. 

^12(^) = ^1=^01, ^2 = 10^, (4.2) 


where g = Q or 
parameter z 


g = A flATil) . 


This i?-matrix is the quasiperiodic function of the spectral 

RUz + l) = Q^^R’l2i^)Qi, 


Ri 2 {^ + r) = exp(—27r^h) A;^ ^ Rui^) 


(4.3) 


on the elliptic curve with periods l,r. At the same time it is the quasiperiodic function of the 
Planck constant h on the smaller elliptic curve with periods 1/N,t/N: 


^fi+r/Ar(^) _ Q^p^_27nz/N) Rizi^) A 2 ■ 


(4.4) 


The property (14.4p can be considered as a consequence of (14.31) and the arguments symmetry 
property 

RUz)=Rf2iNh)Pu, (4.5) 

which follows from flA.12p - flA.13p . In order to verify fl4.4p directly consider 


{T-,K^RUz)(T-,)2 

for some fixed (lA.lTp . Here the lower indices 1,2 are the numbers of the tensor components 
(14.2p . and T~^ = T_.y due to flA.lSp . For the above expression we have: 

^ ^ T h) P—^Pfzt 0 T—aT'y ^ ^ T h.) Pa—'y 0 Ty—a 

Q:€ZiVXZ_/v xWjN 

Now shift the summation indices a ^ a -|- 7. Then from (1A.22P we get 

{T^)y^RUz){T^)2 = exv{27ndrU^)R\t'"\z ). (4.6) 

This answer reproduces fl4.4p for 7 = (1, 0) and 7 = (0,1) (then = Q, = 1/N and = A, 
= t/N respectively). 
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Notice that in the proof of fl4.4l) we used the shift of summation indices (a —)■ a + 7 ). At the 
same time the limits of the sum were not changed because a G xZjv- In fact, we also used 
the periodicity of functions flA.lOp (pa{z,uja + h) in indices considered as discrete variables: 


^ai-\-N,a2 {Z, Ua + UJ{N,0) + h) = ipauci2+N{z, UJa + (^{0,N) + h) = (pa{z, + h) 


(4.7) 


The property (14.7p demonstrates necessity of the exponential factor exp{2mdrUJa) in the dehni- 
tion flA.lOll . 

Symmetric i?- matrix in Mat(A^, ® Mat(ilA C)®^. In this section we assume that N and 
M are coprime integers unless otherwise speciheqj. 


Consider the following expression; 

q;€ Ziv xZm 

where (with the dehnition flA.lOp l 


(4.8) 


^a,a{z, h) = exp ( 2m{z + NoJa)^ + 2mhN^ ) (j){z + NUa , h + OJa) = (4.9) 


( 0^2 \ 

2mhN— j (Pa{z + NUa , h + Ua) 
and for any a = (cxi, 02) G x and a = (di, 0:2) £ x 


UJrv = 


Otl + Oi2'T 

~N 




dl + 0'2'^ 

~M 


(4,10) 


while {Ta} (and {T^}) is the basis in Mat(A^, C) (and Mat(M, C) respectively) dehned as in 
flA.14p - (IA.18p . That is, we mark the Mat(M, C) elements and related quantities (including 

related indices) by tildes. For example, k-^ = exp (^m0ia2 — ^ 2 ^ 1 )/m'^. For coprime N and 

M two sets {f{Noda) , d G Zm x Zm} and {f{ua) , d G Zm x Zm} are equal for a periodic 
function / on the lattice Zm x Zm- 

First, notice that for M = 1 (then d = (0,0)) expression fl4.8|) reproduces the dehnition of 
the Belavin’s Mat(A^, C) i?-matrix fll. 6 |) Ri 2 {z) G Mat(A^, C)®^, and similarly for A^ = 1 we have 
Rf^{h) G Mat(M,C)®T 

The analogues of the quasiperiodic properties (boundary conditions) on the lattice Z © rZ 
are of the form: 

^12,i2(x + 15 ^) = Ql ^ 12 , 12(^5 ^) 5 

(4.11) 

^12 + A ^) = exp(—27r?h) A;^ ^ 12(^1 ^) 


and 


91l2,l2(A^+l) = Qr'^^12,i2(A^) Ql. 

5^12,12(2:, h + r) = exp(- 27 r*z) A^^ ^12,12(2:, K) Af , 


(4.12) 


^The general case is not difficult but requires more additional notations. We consider special case M = N in 
the next section. 
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where Qi = Q In ^ ® 1m, Qi = ^ ® Q ® 1m, and Q, A are the matrices of form 

flA.ldp but of size M x M. 

It is important for us that functions ^a,a{z, K) are periodic in both discrete variables - indices 
a = (ai, 0 ^ 2 ) G Zat x Zat and a = (di, Q! 2 ) G Zm x Zm, i-e. similarly to (14.7p we have 


^q, 2 +A^),q(^, h) fj 2 ),Q(^, h.) 


(4.13) 


and 


^Q,{Ql+M,a2){.^y *1*0 Q2+M) (^, h) *1 *q^(q]^^ 02 ) (^, ^) • (^'l"!) 

These properties allow to shift indices of summations in the same way as it was used in fl4.6p . 
To verify fl4.lip - (14.1211 and fl4.13p - fl4.14p one needs flA.4p . 

In GLat case we had also quasiperiodic behavior f|4.4p on the smaller elliptic curve with 
periods 1/N and t/N. For the symmetric i?-matrix fl4.4p acquires the form: 


^12,i2(^ T ^Af ’ ~ ^12,12(^, Q 2 , 


T 


h 


^i2,i2(^ + N—, h) = exp(-27r?Ar—) A- ^2 


M' 


(4.15) 


and 


^i2,i2(^, ^ T ^) — Qi yii2 i2{z,h)QiQ2, 


(4.16) 


^ 12 , 12 (^ 1 ^+= exp(-27r*—)Ai A- m^2,i2i^,^) M-^2- 

The proof of fl4.15p - fl4.16p is similar to the one given for (14.7p . For example, to get (14.161) consider 


(^7)r^^l2,l2(^,^)(^7)2 = Y1 Y1 ^a.a(^,^)Zo_.^®Va®Zh®T_o. 

(a€ Zjv xZ/y/- aG 

Then, shift the summation index as a —)■ a -|- 7 and notice that 

*l’o+7,o(^, h) = exp{2mzdrUJ^)k^^^^a,a{^, h + uj^). 

Plugging 7 = (1, 0) (or 7 = (0,1)) into the obtained expression one gets the upper (or the lower) 
line of fl4.16p . 

Quadratic relations. Let us prove that the symmetric i?-matrix obeys the relations similar 
to dni), (\n^ and (Op. 


Proposition 4.1 The symmetric R-matrix (4:_8) satisfies the following set of properties: 

1. Skew-symmetry: 

^ 21 , 2 i(~^, ~^) — ~^ 12 ,l 2 (^, • (4-17) 


2. Associative Yang-Baxter eguation: 


^12,12 ^23,32 — ^13,32 ^12,13 + ^23,31 ^13,12 


(4.18) 


where = 91^6 ,Ha - hi). 
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3. Unitarity: 


^12,12 ^21,12 ~ ® ® IjU ® 1 m N M [p[Nh) — p{Mz)) . 


(4.19) 


or 


^ 12,12 ^ 12,21 ~ Iat 0 Iat <8) 1m 0 1m N M {p{Mz) — p{Nh)) . 


(4.20) 


Proof: The first property simply follows from the dehnition fl4.8p - fl4.9p : 

OiG Ziv xZ_/v oE Zm xZm 

Changing snmmation indices as a —)■ —a, a —a and nsing $_« _«(— 2 ;, —h) = —^a,a{z, h) we 
get (14.171) . 

The second property follows from (IA.19P (for both k and k) and the Fay identity flA.SD 
written in terms fnnctions (14.9p as 


(4.21) 


<l'„,a( 2 ;, h) = 

= ®/3,a+/3(^ + 3 ) ®a-/3,a(^, h - T]) + T] - k) + W, K) 

with z = Zi — Z 2 , w = Z 2 — Z 3 , h = hi — ^2 and t] = h^ — h^- 
To prove the nnitarity (I4.19P let ns write its l.h.s. explicitly 

^12,l2(^) ^21,12(~^) ^a,a{z, h) $^^^( — 2:, H) 0 T^-q 0 0 T_-_^ 

a,/3,a,/3 

and snbdivide this snm into fonr parts as 

E y E y E y E , 

a,/3,a,/3 a=p ,a=—^ a=^ a^/3 


The hrst snm reprodnces the answer dne to (IA.7j) and (lA.l Hfl: 


^ <h«,a(^, h) (-^, h) = '^'^p{h + Ua) - p(z + Nua) = N‘^M‘^{p{Nh) - p{Mz)). 

a=l3,a=-y “ “ 

Each of three other snms eqnals zero. The vanishing of the second and the third snm is proved 
in the same way as it is made for nnitarity of R^{z) in the Appendix. Consider the fonrth snm. 
Denote 'y = a — (3,^ = a + (3. Then, nsing (14.211) 


a^/3,aj^-y 77 ^ 0 , 7^0 


X 


^ 1 Y1 ^07 *^/3,7(0, K) <F^,a(z, 0) + ^ <^’^,7(0, h) 0) 

a,/3 


^We also use here that ^ p{z + NPa)=J2 p(z + oJa) for coprime N and M. 
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The latter two terms cancel each other after changing notation a := /3, := —a in the last 

term. This hnishes the proof of fl4.19p . At last, fl4.20p follows from fl4.19p and fl4.17p . ■ 

Cubic relations. The next step is to derive cnbic relations of (12.ip - (12.41) type. Recall that (12.31) 
is the most general in that list. In a similar way one can obtain its analogne for the symmetric 
R-matrix fl4.8p . First, mnltiply fl4.18l) by 9323,13 from the left: 

fr^23,i3 fr^l2,i2 fr^23,32 = ^^23,13 ^^13,32 ^^12,13 + ^^23,13 ^^23,31 ^^13,12 • 

Second, change indices 2 o 3 in fl4.18p and mnltiply it by 933332 from the right: 

fr^l3,i2 fr^32,32 fr^23,32 = ^^12,32 ^^13,13 ^^23,32 + ^^32,31 ^^12,12 ^^23,32 • 

Combining these two eqnations and nsing the properties of skew-symmetry fl4.17l) and nnitarity 
flCTll - frai we get 

fr^l2,32 fr^l3,i3 fr^23,32 ~ ^^23,13 ^^13,32 ^^12,13 + A^ Af 93]^3,i2 (^P(A^% 2 ) ~ ) (4.22) 

where h-^ = ha — fi^- When M = 1 (I4.22p coincides with (12.31) . where h = hi —and r] = hg — 

In particnlar case of hg = {hi + h^)/^ we get Yang-Baxter like eqnation 

fr^i2, 32 (^ 12 ) 93^3 ^ 3 (^ 13 , h) 9323 32 ( 2 ^ 23 , ~ fr^23, 13 (^ 23 , h) 9333 32 (^ 13 , h) 933243 ( 2 ^ 12 , h) (4.23) 

with h = hjg = h^ 2 - 

Another particnlar case hi = h^ in fl4.22l) leads to analogne of fl2.2l) : 


^^21,23(^21, h) 9333 33(^33, h) 9332 23(^32, h) -|- 9323 33(^23, ^31,23(^31, ^) 9332^33(^12, h) — 

= -u ® Ijv ® hs N^M^p'{Nh ). 


( 4 . 24 ) 


A similar derivation leads to the cnbic eqnation of ( 14 . 221 ) type with h and z variables inter¬ 
changed: 

fr^ 32 ,i 2 fr^l3,i3 fr^ 32,23 ~ ^^ 13,23 ^^ 32,13 ^^ 13,12 + N M 933243 ^p(M 2 ; 23 ) ~ ^(^^2:13)^ , ( 4 . 25 ) 

where Zab = Za - Zb- 

Now we can easily obtain the cnbic Yang-Baxter like eqnations of ( 12 . 41 ) type. It is based on 
the skew-symmetry of ^(Yhgg) — p{Nhi^ with respect to interchanging of the argnments only. 
Therefore, from fl 4 . 22 l) we have: 


fr^i2,32(^2, ^32) 9333 33(^13, hgg) 9323 32(^23, ^32) + fr^i2,32(^2, ^13) 9333 33(^13, ^32) 9323 32(^23, ^13) 


~ fr^23,13(^23, ^i3)fr^l3,32(^3, ^32)fr^l2,13(^2, ^13) + 9323 13(^23, ^32)fr^l3,32(^3, ^13)9332^13(2:12, h.32) 
In the same manner from ( 14 . 251 ) we get: 

fr^32,12(^32, ^12) 9333 33(2:13, hjg) 9332 23(2:32, hgg) -|- 933212(2:13, ^32) 9332 33(2:13, hjg) 9332 23(2:13, ^23) 


~ fr^i3,23(^13, ^23)^^32,13(^32, ^13)^^13,12(^13, ^12) + ^^13,23(^32, ^23)^32,13(^13, ^13)9313^12(^32, ^12) 
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Relation to Baxter-Belavin GLjvm -R-matrix. Multiply the symmetric i?-matrix by NP 12 ® 

Im <8 Im- 

q:,/5e ^iv xZiv aEZjpfXZjpf 


IXT 31 


X ^^0,7 QW{NhdrUJa)^a{,Z + h + CJo) ® T_^ 0 Tq, ® T_a G 


a, 7 ,a 


X o;^ + (ha + 0 T_^ ®Ta® T_a ■ 


7,0; 


Since M and N are coprime the summation over oj^ + Ua is equivalent to the summation over Ua, 
where a G Zmn x Zmn- Therefore, we get the expression which is the GLtvm Baxter-Belavin’s 
R-matrix written in special basis. 


Particular case M = N. In this particular case the function ^a,a{z^h) fl4.9p is simplihed. 
Using flA.4p it is easy to get that 


^a,a{z, h) = exp( 27 r? 


0102 — O 2 O 1 

N 


)ipa{z,u:a + h) = Kl^^ipa{z,Ua + H) 


5 Appendix 

5.1 Elliptic functions 

We deal with the following elliptic functions: 


(l){z,u) = 


'&'{0)'&{z + u) 
'd{z){}{u) 


EAz) = 


'd'{z) 


^{z) ’ 

where i9(z) is the Riemann theta-function 


E 2 (z) = -d^Ei(z) = p(z) 


ir'(o) 

3 'd'(O) ’ 


'd(^) = d{z\T) = X^^P + 27r*(2; -h ^)(/c -h j 

fcez ' 

which has simple zero at 2 ; = 0, and p(z) in flA.2p is the Weierstrass p-function. 
The function <^(x, y) = x) has the quasiperiodic properties 

(l){x + 1,2/) = (l){x, y ), (j){x + T,y) = exp{-27ny)(j){x, y) 


and satishes the Fay type identities: 

0 (x, n) 0 ( 2 /, w) = 0 (a; — //, u)(l){y, u + w) + (j){y — x, w)(l){x, u + w), 

(p{x, z)(p{x,w) = (p{x, z + w){Ei{x) + Ei{z) + Ei{w) - Ei{x + + w )), 

(j){z, u)(j){z, -u) = p{z) - p{u) = E 2 {z) - E 2 {u) . 


(A.l) 

(A. 2 ) 


(A.3) 


(A.4) 
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(A.5) 

(A. 6 ) 

(A.7) 















These relations lead to 


(p{z^ Ui — v)(j){w, U2 + — W,v) — (j){z, U2 + v)(f){w, Ui — v)(l){z — W,Ui — U2 — v) 

= (j){z,Ui)(j){w,U2){Ei{v) - Ei{ui - U2- v) + Ei{ui - v) - Ei{u2 + v)) 


(A.8) 


and 


(j){z, U — v)(f){w, v)(j){z — W,v) — v)(l){w, u — v)(j){z — w,u — v) 


= 4 ’{z,u){p{v) - p{u - v)). 
We also need the following set of fnnctions: 


Paiz, Ua + E) = exp (^2tii^z'^ (f){z, Ua + K), a G 


(A.9) 


(A.IO) 


These fnnctions satisfy the set of identities flA.5P - flA.9p because the exponential factors are 
canceled. 


Averaging p-function: 


+ h) = N‘^p{Nh) 

xIa ]^ 


(A.ll) 


The following formulae (of hnite Fourier transformation type) are also useful in derivation 
of i?-matrix identities (see d]): 

1 ^ 

^ + ]^) = > VyGZ''^ ('A.12) 

XZjv 

or (exchanging the arguments) 

1 ^ 

Y1 ^l,'y¥’a{z,Ua, + h) = ip^{Nh,u;^ + —), V7GZ^2^ ('A.ia) 

xZjv 

where is from flA.lSp . 


5.2 Baxter-Belavin i?-matrix 

Finite-dimensional representation of the Heisenberg group. Consider the following 

N X N matrices: 

Stt'Z 

Q,A G Mat(W,C) : Qki = hiexp{—k), Eki = 6 k-i+i= 0 modN , fc, / = 1,..., iV, (A.14) 

= A^ = 1 , (A.15) 

exp(27r*^^) Q'^iA '^2 = , 71 , 72 G Z. (A.16) 

Introduce 

T, := T,,,, = exp( 7 r*^) , 7 = ( 71 , 72 ) G Z"^ . (a.17) 

The subset of with 71,72 = 0,..., W — 1 is the basis in Mat(W, C). From flA.lbp we have 

TaTp = Ka,pTa+p , = GXp (A «2 “ l^ 20 tl)j , (A.18) 
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where a + (3 = (cti + /3i, 02 + (32)- Obviously, 


l^a,l3^P,a 1 ) q:,/3 1^/3,a i ^a,/? ^Q!+/3,/3 l^a,)3+a • 


Notice also that 


5^<, = ]VXo. 


(A.19) 

(A.20) 


where 0 = (0, 0). 

For N = 2 the basis (IA.17P in 2 x 2 matrices coincides (up to signs) with the Pauli matrices 
basis cTj endowed with uo = 12x2- In this case flA.lSD can be re-written as crjak = o’Qdjk + iSjknO'n 
together with cTjCTo = a^aj = aj and ag = ao- 

i?-matrix is dehned a^ 


+ h) Tq, 0 T_a G Mat(iV, C)®^ , 


where 


ZjY X L 


0^2 0^2 0^2^ 
^Paiu.oJa + h) = exp( 27 r*—u) h), a;« =- — - 


(A.21) 


(A.22) 


Writing R^bi^) we mean that Ta®T_a in 1IA.21I1 is replaces by 10...10Ta 01...1 0 T_q, 0 1... 01 
in such a way that and T_q, acts in the a-th and 6-th components of Mat(iV, If the 

number of components equals n then G Mat(iV”',C). With the above dehnition we get the 
correct relation 

Rla{z) = Pah Rlb{-^) Pah , (A.23) 

where P 12 is the permutation operator in Mat(A^, C)®^. Notice also that P12 has the following 
form in the basis (1A.17P : 

= NPu, 


The classical limit h —)■ 0 is dehned as 

Rui^) = ® 1 + ri2{z) + hmi2{z) + 0{h^). 

Then the quantum Yang-Baxter equation (ll.3p provides the classical one 
[rah, Tac] + Vac, Tbc] + Vab, The] = 0 , = TabiZa “ Zb) • 

Condition (ll.7p leads to 


^ah '^ha : ’^ab ^6a : 

while (ll.Sp gives expresses niab in terms of Vab- 

2 mi2{z) = rl^iz) - 1 0 lN‘^p{z). 

The answer is similar to h} term in the local expansion of the Kronecker function near h = 0: 

z) = tr^ + E^{z) + h {El{z) - p{z))/2 + 0{te). 


(A.24) 

(A.25) 

(A.26) 

(A.27) 

(A.28) 


®Let us mention that original definition differs from (IA.21I) . It was defined in terms of theta functions with 
characteristics 1/A^ and normalized as R 12 R 21 = 1 ® 1- Its representation in the standard basis of Mat(7V, C) 
{(^b)fei = ^^kSji} was suggested in [12]. 
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(A.29) 


At the same time the Fay identity for /2-matrix fll.2p gives 

Tab Tac “ Uc Tab + Tac He = TUab + rUbc + mac ■ 

Plugging niab from flA.28p into flA.29p one gets the non-abelian analogue of 

(^Eii^Za Zb) E-i(^Zb Zc) E\{^Zc ^a)) ^{,^a ^b) ^{,^b ^c) p(^c ^a)- 

In the end let us prove the unitarity property (ll.Sp (this proof was skipped in [I3l E]). 
Unitarity. Proof of fll.81) . It is convenient to rewrite fll.8p using fll.7p as 

RUz)RT 2 \z) = -N\p{Nh) - p{z)) 10 1. 

The proof is achieved by direct calculation: 

R^2iz)Ri2iz) = ® T_aT_p(pa,{z,UJa + h)(pf}{z, UJfj “ 

a,13 a+P =0 a+13^0 

The hrst term reproduces the answer due to flA.7p and J2a p{^a + h) = N‘^p{Nh). So we need 
to prove that the second term (where a + (3 ^ R) equals 0: 

^ ^ T^0T_.,<^(^^(;2,n;.,)(/?i(;2) + /?i(a;„-fh) + /?i(n;^-h)-/?i(^-l-a;.,)) 

a+Z^T^O 77^0 a+/3=7 

The hrst and the last terms in the brackets are independent of a, (3. Each of them gives 0 since 
i^a ,(3 = i^a,a +/3 = i^a,^ and but the sum does not contain 7 = 0 term. Finally, 

we are left with 

-t- h) -|- Eiiojp — h)) = K^ .y{Ei{uJa + h) + Ei{ujj—UJa — R)) =0 

Q;+/3=7 Oi 

because the summation index of the second term can be shifted as a ^ a -|- 7, and due to 
Ei{x) = -Ei{-x). 

More details and properties can be found in [191 [12] and [Tl| . 


5.3 GLyr Sklyanin algebra 


The Sklyanin algebra [21] has extension to GLat. One possible approach is given in [T^. The 
most natural way is to use the exchange relations (as in |21jh 


Proposition 5.1 Consider the Lax operator / O.ill) and R-matrix U.6]) . The quantum exchange 
relations 

R\ 2 {z - w)L\{z)Ll{w) = LI{w)L\{z)RUz - w) 
are equivalent to x relations numbered by a, j3 E (for the component Ta 0 Ty): 


E 


,a—l3Sa—'yS(3+^ 1q;,/3,7 


0 , 


where the structure constants are given by 

for (3^0: fa, ^,7 = El ( 1^7 + h) — Ei + h) + Ey {uJa-^ + h) — El (ci ;^+7 + h) , 

for /3 = 0 : f^ = p{uj^ + h) - p(ua-y + h) . 


(A.30) 


(A.31) 
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For N = 2 relations flA.30|) reproduce GL2 Sklyanin algebra in its original form [2l]. Indeed, 
for a, 7^ 0 and a ^ (3 the structure constants are related: f^ ^ p ~ ~^a 0 0 

fa,/3,a = fa,/3,o- Therefore, 

[Sa, 5/3]- = KaASa+0, 5o]+ . (A.32) 

Here we also used that for a,[3 and a ^ (3 Ka,0 = ±*, i.e. Ka,0 = —1^0,a- To get the second 
commutation relation introduce 


and notice that f^ o,a ~ 


p1 = p{h + ujA - p{h) 

_eh ch _ _fh Thpn 

fo,0,0; ‘■a,0,0 ~ ‘■a,O,a-0- nen 


pI [^0, A+0]- = -Ka,0{pa - pj)[5a, A]. 


(A.33) 


(A.34) 


for some (3 ^ a,0. Notice also that the same commutation relation appears from a = 0, /5 7^ 0 
component of flA.30p but in a different form: 


|So, S^l+ , 

where = Ei{h + Ua) — Ei{h) — Eiipja)- The latter is the same as flA.34D due to 

Kg- ^ Pa - P0 

^a+0 Pa+0 

in GL2 case. 
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